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atoms and show that it is possible to obtain a large degree of steady-state entanglement 
in the system. Necessary conditions for the entanglement are nonclassical two-photon 
correlations and nonzero collective decay. It is shown that nonidentical atoms are a 
better source of stationary entanglement than identical atoms. We discuss the optimal 
physical conditions for creating entanglement in the system, in particular, it is shown 
that there is an optimal and rather small value of the mean photon number required for 
creating entanglement. 
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1. Introduction 

Entanglement between separate quantum systems is one of tlie key problems in quantum 
mechanics. A number of interesting concepts and methods for creating entanglement 
have been proposed involving trapped and cooled ions or neutral atoms mEllSlllllSllElIZllH! • 
Of particular interest is generation of entangled states in two-atom systems, since they 
can represent two qubits, the building blocks of the quantum gates that are essential to 
implement quantum protocols in quantum information processing. It has been shown that 
entangled states in a two-atom system can be created by a continuous driving of the atoms 
with a coherent or chaotic thermal field |H|FHll()piipi2j . or by a pulse excitation followed by 
a continuous observation of radiative decay P^14pi5j . Moreover, the effect of spontaneous 
emission on initially prepared entangled state has also been discussed [THlll7pi8|IT^ . These 
studies, however, have been limited to the small sample (Dicke) model [201 or the situation 
involving noninteracting atoms strongly coupled to a cavity mode. The difficulty of the 
Dicke model is that it does not include the dipole-dipole interaction among the atoms 
and does not correspond to realistic experimental situations of atoms located (trapped) 
at different positions. In fact, the model corresponds to a very specific geometrical 
configuration of the atoms confined to a volume much smaller compared with the atomic 
resonant wavelength (the small-sample model). The present atom trapping and cooling 
techniques can trap two atoms at distances of order of a resonant wavelength pi j l ^ l^ , 
which makes questionable the applicability of the Dicke model to physical systems. 

Recently, we have shown j21] that spontaneous emission from two spatially separated 
atoms can lead to a transient entanglement of initially unentangled atoms. This result 
contrasts with the Dicke model where spontaneous emission cannot produce entanglement 
from initially unentangled atoms fOl^]- We have also found |2S1 analytical results for 
two measures of entanglement and the relation between them for the two-atom system 
radiating by spontaneous emission for quite broad range of initial conditions. 

In this paper we study the creation of a stationary entanglement in a system of 
two identical as well as nonidentical two-level atoms separated by an arbitrary distance 
ri2 and interacting with a squeezed vacuum. The squeezed vacuum appears here as a 
source of nonclassical two-photon coherences, essential for the creation of the stationary 
entanglement. We use the master equation to describe the evolution of the system and 
find the steady-state solutions for the atomic variables. We present analytical results 
for concurrence which is well known and calculable measure of entanglement. We find a 
surprising result that non-identical atoms with significantly different transition frequencies 
can exhibit a larger entanglement than identical atoms. Under some conditions, the 
nonidentical atoms can be maximally entangled with the value of the concurrence equal 
to unity. 
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2. Master equation 

We consider a system of two two-level atoms at fixed positions ri and r2 and coupled 
to the radiation field, whose the modes are in a squeezed vacuum state. Each atom has 
energy levels \gi) and |ej) (i = 1,2) such that -Eg, — Eg. = huji, transition dipole moment 
/Xj, which we assume equal for both atoms. 

We analyse separately the dynamics of identical and non-identical atoms. In the case 
of nonidentical atoms, we assume different transition frequencies ui and uj2 such that 
A = {lj2 — t^i)/2 <^ ti^o = ("^1 + '^2)72, so that the rotating-wave approximation can be 
applied to calculate the dynamics of the system. 

The system can be described by the reduced density operator p which, in a 
Schrodinger picture, satisfies the master equation j^EI 

§^ = - ^ E (1 + N) {pStSr + S^Srp - 2SrpSt) 

1 2 

r.,iV {pSrS; + SrS]-p - 2S^pSr) 



1 2 

+ ^^JM* {pSrS. + SrSrp - 2SrpS.) e 



2 2 

[s^, p] - ^ E [stsi,p] . (1) 

i=l i^j 

Here, and S~ are the raising and lowering operators, respectively, of the ith atom, 
and M = \M\ exp (i0s) characterise squeezing such that |Mp < A^(A^ + 1), where the 
equality holds for a minimum-uncertainty squeezed state, (ps is the squeezing phase and 
ujs is the carrier frequency of the squeezed vacuum. 

The parameters Tij, which appear in equation ((T)), are spontaneous emission rates, 
such that 

is the spontaneous emission rate of the ith atom, assumed to be equal for both atoms, 
and 

21 sin (Via) 



ri2 = r2i = ^Vf|[l-(M-ri2)'] 



+ [1 - 3 (A ■ ri2)'] 



kori2 

cos (fcori2) _ sin (fcori2) 

. (Vl2)^ (Vl2)^ 



(3) 
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are collective spontaneous emission rates arising from the coupling between the atoms 
through the vacuum field [SHlEniEZI, and 

3 r ri - n2i cos(A;ori2) 



f^i2 = ^r|-[l-(A-ri2f] 



+ [1 - 3 (A • r.,)'] 



sin(A;o?"i2) , cos (/cori2 



+ 



represents the vacuum induced coherent (dipole-dipole) interaction between the atoms. 

In the expressions Q and (jH), fi and ri2 are unit vectors along the atomic transition 
dipole moments and the vector = r2 — ri, respectively, and ko = uq/c. Later on, we 
will assume that the atomic dipole moments fx are perpendicular to the vector joining 
the two atoms. 

The collective parameters ri2 and Q12, which arise from the mutual interaction 
between the atoms, significantly modify the master equation of a two-atom system. The 
parameter ri2 introduces a coupling between the atoms through the vacuum field that the 
spontaneous emission from one of the atoms infiuences the spontaneous emission from the 
other. The dipole-dipole interaction term Q12 introduces a coherent coupling between the 
atoms. Owing to the dipole-dipole interaction, the population is coherently transferred 
back and forth from one atom to the other. 

The two-atom system can be described in the basis of product states of the individual 
atoms 

|l) = \9i) ® 1^72) = \g), 

|2) = |ei) (g) |e2) = |e) , 
|3) = \gi) ® lea), 

|4) = |ei) ® \g2) , (5) 

where \gi) and \ei) (for i = 1,2) are the ground and excited states of the individual atoms. 
In this basis, the two-atom system behaves as a single four-level system whose the density 
matrix can be written as a 4x4 matrix. Due to the presence of the dipole-dipole interaction 
Q12 it is often convenient to introduce collective atomic states that are eigenstates of 
the system of two identical atoms including the dipole-dipole interaction [201 ■ They are 
symmetric and antisymmetric superpositions of the product atomic states |3) and |4), 
given by 

\s)= ^(|3) + |4)), 

l«)= -^(|3)-|4)), (6) 

and the states \g) and |e) remain unchanged. 

On introducing the collective states © and using the master equation we are able 
to write down equations of motion for the matrix elements of the density matrix for two 



(4) 
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atoms in a squeezed vacuum. We treat separately two cases of identical and nonidentical 
atoms. 

2.1. Identical atoms 

For identical atoms separated by an arbitrary distance ri2 and interacting with a squeezed 
vacuum field of the carrier frequency Ug = ujq, we obtain the following set of coupled 
equations of motion [2E] 

Pee = - 2r (AT + 1) Pee + AT [(F + Tu) Pss + (F - F12) Paa] + Fi2|M|p, , 
Pss = (F + F12) {N - (3iV + 1) p,, - Np,a + Pee ' |M|pJ , 
Paa = (F - F12) {N - (3iV + 1) Paa - Np,, + p,, + |M|pJ , 

p, = 2F12IMI - {2N + 1) Fp„ - 2\M\ [iV + 2F12) Pss - (r - 2F12) Paa] , (7) 

where p„ = Pe3exp(— z^^) + pgeexp{i(j)s). It is seen from equation ((7j) that the evolution 
of the populations depends on the two-photon coherences p^g and Pge, which can modify 
population distribution between the collective states. The coherences can also create 
superposition (entangled) states involving only the ground \g) and the upper |e) states. 
The evolution of the populations depends on F12, but is completely independent of the 
dipole-dipole interaction ^12. 

The steady state solutions of equations ((Zj) depend on whether F12 = F or F12 7^ F. 
For two atoms separated by an arbitrary distance ri2, F12 7^ F, and then the steady-state 
solutions of equations ((Tj) are 

_ [{2N + If - 4|Mp] + |M| 

~ {2N + 1)4 - 4|M|2 [{2N + 1)2 - 723] ' 

_ N{N + 1) [{2N + 1)2 - 4|M|2] + |M|27i2(7i2 - 2) 

~ (2Ar + l)4-4|M|2[(2iV+ 1)2-72^] 

^ N{N + 1) [(2iV + 1)2 - 4|M|2] + |M|27i2(7i2 + 2) 
(2Ar+l)4-4|M|2[(2Ar+l)2-722] 

2(2iV+l)|M|7i2 . . 

(2Ar+l)4-4|M|2[(2Ar + l)2-722] ' 

where 712 = F12/F is the dimensionless collective damping parameter. This result 
shows that all the collective states are populated in the steady-state even for small 
interatomic separations (712 ~ 1). For large interatomic separations 712 ~ 0, and then 
the symmetric and antisymmetric states are equally populated. When the interatomic 
separation decreases, the population of the state \a) increases, whereas the population 
of the state |s) decreases and pss = for very small interatomic separations. This 
effect results from the enhanced (F + F12) damping rate of the symmetric state, as it 
is seen from equation ((Tj). The two-photon coherences, represented by p^, affect the 
population distribution only when both \M\ and 712 are nonzero. The coherences are 
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crucial for getting entanglement in the system. Of particular interest are population 
distributions for maximally squeezed fields with \M\ = ^yN(N + 1). In this case, the 
factor (2A^ + 1)^ — 4|Mp = 1, and then the solutions (jH)) take a very simple form 

_ + N{N + 1) 

iV(iV+l)(l- 7,2)2 



Pa 



l + 4iV(iV + l)(l + 722) ' 

iV(iV + !)(! + 7,2)^ 
l + AN{N + l){l + ^f,) ' 



Pu — 1 , ^ AT/ AT , , 2 \ ' y ) 



2y/N{N + l){2N + l)^u 
l + AN{N + l)il + jj,) 

The solutions (jHl) and will be used for calculation of the degree of entanglement present 
in the system. For further reference it is important to note that the sum of the populations 
Pss + Paa tends to 0.5 as the quantity A^(A^ + 1)(1 + 7^2) becomes much greater than one, 
which means that for large values of the mean number of photons one half of the 
population goes eventually to the states |s) and \a). 

2.2. Nonidentical atoms 

The population distribution is quite different when the atoms are nonidentical with 
A = {ijj2 — ijJ\)l2 7^ 0. As before for the identical atoms, we use the master equation 
and find four coupled differential equations for the density matrix elements with time- 
dependent coefficients oscillating at frequencies exp(±zAt) and exp[±2i(a;s — ti;o)t-l-0s]. If 
we tune the squeezed vacuum field to the middle of the frequency difference between the 
atomic frequencies, i.e., ujg = {ui+uj2)/2, the terms proportional to exp[±2i(a;s— ci;o)t-|-0s] 
become stationary in time. None of the other time dependent components is resonant with 
the frequency of the squeezed vacuum field. Consequently, for A ^ F, the time-dependent 
components oscillate rapidly in time and average to zero over long times. Therefore, we 
can make a secular approximation in which we ignore the rapidly oscillating terms and 
obtain the following equations of motion 

Pee = - 2F (iV + 1) Pee + A^F (p,, + p,,) + Fi2|M|p„ , 
Pss =T[N-{3N + 1) Pss - Npaa + Pee] - Ti2\M\pu , 
Paa = T[N-{3N + 1) Paa ' Np^s + Pee] ' Ti2\M\p^ , 

p, = 2F12IMI - {2N + 1) Fp„ - 4Fi2|M| {pss + Paa) . (10) 

The steady-state solutions of equations (fTO|) are 
1 f (2A^- 1) 1 



Pe 



4 1 2N+1 [{2N + 1)' - 4|M|2 7^2] 
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p = 21^1^^^ (11) 

(2iV+l) [(2iV+l)'-4|M|2 7y ■ 

Equations are quite different from equations (jH)) and show that in the case of non- 
identical atoms the symmetric and antisymmetric states are equally populated. This 
fact will be crucial in the entanglement creation in the system and results from the equal 
damping rates of the symmetric and antisymmetric states, as it is seen from equation (|T0|) . 
For maximally squeezed vacuum with |M| = ^^/W(7^^^^1) the solutions (fTT|) simplify to 

1 (2N -1 1 

Pee = T i TTTT— r + 



4\2A^+1 1 + 4iV(iV + 1)(1 -71^2) 
iV(iV +l)(l-7?2) 



p. 



2v/iV(iV + 1)712 

(2iV+l)[l + 4iV(iV + l)(l-7?2)] ■ ^ ' 

From equations (fT^ it is evident that for the Dicke model, for which 712 = 1, the 
populations of the symmetric and asymmetric states both are zero and p„ tends to unity 
for large A^. In real situations the separation of the atoms is nonzero, and we have always 
7i2 < 1, which means that for N{N + 1)(1 — 7^2) ^ 1 the populations Pss and Paa both 
approach value 0.25, i.e., pss + Paa ~ 0.5. That is, for sufficiently large N one half of 
the population is transferred to the block spanned by the symmetric and antisymmetric 
states, similarly to the identical atoms. There is, however, one essential difference between 
the identical and nonidentical atoms, which is the scale for the saturation - much slower 
for nonidentical atoms. For 712 very close to unity, very large intensities N are required 
to have A^(A^ + 1)(1 — 7^2) ^ 1, and consequently pss + paa ~ 0.5. Thus, for small N and 
7i2 7^ 1, the populations pss and paa are very small, pss = Paa ~ 0. This fact will have 
important effect on the entanglement creation in the system of nonidentical atoms. 



3. Steady state entanglement 

To assess how much entanglement is stored in a given quantum system it is essential to 
have appropriate measures of entanglement. A number of measures have been proposed, 
which include entanglement of formation j2Zj, entanglement of distillation [2H], relative 
entropy of entanglement j^H] and negativity IHOIISIIESESI- For pure states, the Bell 
states represent maximally entangled states, but for mixed states represented by a 
density matrix there are some difficulties with ordering the states according to various 
entanglement measures; different entanglement measures can give different orderings of 
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pairs of mixed states and there is a problem of the definition of the maximally entangled 
mixed state jSUES]. 

Here we use the concurrence to describe the amount of entanglement created in a two- 
atom system by the interaction with the squeezed vacuum. The concurrence introduced 
by Wootters [2Z| is defined as 

(13) 



C = max (^0, V Ai — V A2 — V ^3 ~ V 
where {Aj} are the the eigenvalues of the matrix 

R = pp (14) 
with p given by 

p = (T^ (g) p* ^ ax , (15) 

ax is the Pauli matrix, and p is the density matrix representing the quantum state. The 
range of concurrence is from to 1. For unentangled atoms C = whereas C = 1 for the 
maximally entangled atoms. 

In the basis © of the product atomic states the density matrix for two atoms in the 
squeezed vacuum has in the steady state the following block form 



(16) 



has the form 





Pii 


P12 





\ 




P21 


P22 
















P33 


P34 


v 








P43 


P44 / 


Tr p 


= 1. 
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matrix 
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P22 


P12 
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P21 


pll 
















P44 


P34 


\ 








P43 


P33 / 



(17) 



and the square roots of the eigenvalues of the matrix R given by ()14|1 are the following 

{v^}= {VP11P22 - IP12I, VP11P22 + IP12I, VP33P44 - IP34I, A/P33P44 + |P34| } • (18) 

Depending on the particular values of the matrix elements there are two possibilities for 
the largest eigenvalue, either the second term or the fourth term in (fTHj) . The concurrence 
is thus given by 

C = max{0, Ci, C2}, (19) 

with 



Cl = 2 (IP12I - VP33P44 ' 
C2 = 2 (|p34| - VP11P22 . 



(20) 
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and we have two alternative expressions for the concurrence depending on which of them 
is positive. 

In terms of the collective atomic states \g), |e), |s) and |a), the expressions for the 
concurrence l^20\i take the form 

Ci = 2 Ipgel - V {pss + PssY - ipsa + Pasf ) , 

C2 = \Pss - Paa + Psa - Pas | - 2 y/PggPee ) ■ (21) 

Having the concurrence expressed in terms of the density matrix elements, we can apply 
the steady-state solutions (jH)) or |TT| and obtain analytical results for the stationary 
concurrence. We will discuss the results separately for identical and nonidentical atoms. 



3.1. Identical atoms 

In case of identical atoms the steady-state solutions (jH]) are valid, and we find 

Cl = \pu\ - (Pss + Paa) 

_ (2iV + 1)|M| 712 - |Mp7f, - NjN + 1) [{2N + 1)^ - 4|Mp] 

(2Ar + 1)4 -4|M|2[(2Ar+ 1)2-722] • ^ ^ 

It turns out that C2 is always negative, so the only contribution to the concurrence comes 
from Cl, and then C = max(0,Ci). It is clear from (f^^ that the perfect entanglement in 
the system, i.e., the value of concurrence equal to unity, can be achieved when = 1 
and Pss + Paa = 0. This is generally impossible for identical atoms because there is 
always some population stored in the states \s) and \a). However, it is possible to obtain 
some degree of entanglement in the system for appropriately chosen values of ri2, and 
|M|. The concurrence, measuring the degree of entanglement, depends on the interatomic 
distance through the collective damping parameter 712, and the degree of the two-photon 
coherences |M|. It is easily to show from equation ()22j) that there is no entanglement 
possible for \M\ < N, i.e., for classically correlated fields. For a quantum squeezed field 
with maximum correlations \M\ = ^/iV^iVH-Ty, the concurrence can be written as 



r -g / .r(.r , (2iV + 1) 712 - ViV(iV + 1) (1 + 71^2) .^S^ 

Cl - 2 VNiN + 1) TTim^WT^) ■ ^ ^ 

In this case, Ci can be positive. To show this, we plot in Fig. ^ the concurrence for 
two identical atoms in the maximally squeezed vacuum as a function of the interatomic 
distance ru and the mean number of photons A^. It is evident from Fig. ^ that there is a 
range of values of ^12/ A for which C is positive. The maximum of concurrence is obtained 
for ri2 ~ when the atoms are very close to each other. The values of concurrence 
decrease as the interatomic distance increases and reduces to zero at ~ A/2, but we 
can observe revival of concurrence for longer interatomic distances, although the next 
maximum is much weaker. It is interesting that the maximum of concurrence appears for 
not very high values of the mean number of photons < 0.1. It is easy to check, from 
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Figure 1. Concurrence C for two identical atoms as a function of the interatomic 
distance ri2/A and the mean number of photons A'' for |M| = ^/N{N + 1) 



that for = 1 the number of photons Nj^ax for which Ci reaches its maximum is equal 

to N^ax — i^J (1 + a/2)/2 — 1)/2 Ki 0.049. The fact that the maximum of concurrence 
appears for moderate values of the mean photon numbers can be important from the 
experimental point of view as the present sources of squeezed fields can produce quantum 
squeezed fields of intensities < 1. 

The steady-state entanglement and its presence for only quantum squeezed fields of 
small intensities N is associated with nonclassical two-photon correlations characteristic 
of the squeezed vacuum field. To show this, we introduce a parameter 
IMI 

rC=y, (24) 

which characterises two-photon correlations normalised to the intensity of the squeezed 
field. For classical fields, \M\ < N, and then TC < 1 for all N. For a quantum squeezed 
field with |M| = y/N{N + 1), the parameter becomes 

TC = fT^, (25) 
which is always greater than one. The result is a strong two-photon correlation, which 
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is greatest for < 1. Thus, the non-classical two-photon correlations are significant for 
< 1 and lead to a large entanglement in the system. 

3.2. Nonidentical atoms 

In the case of nonidentical atoms with A 7^ 0, the steady-state values for the 
density matrix elements are given in equation |TT| . As above for identical atoms, the 
concurrence ()19|1 can be expressed by the formula piU) which, with the solutions (fTT|) . 
leads to 

2|M||7i2| 



Cl = \Pu\- (pss + pa 



i2N+l) [(2Ar+lf -4|MP7y 

-111 ^ 1 (26) 

2 I (2Ar + l)2-4|M|2 72j 

Similarly to the case of identical atoms, C2 is always negative. Moreover, Ci is always 
negative for |M| < A^ independent of 712. Thus, entanglement is possible only for quantum 
squeezed fields which for the maximum correlations |M| = a/ A^(A^ + 1) gives 



. ^,,-{2N + l)^N{N + l){l--ff,) 
i2N + l)[l + AN{N + l)il-^f,)] 



C, = 2 VN{N + 1) '^';J^ /LrZ ~2 ^? • (27) 



Equation (jTTj) is significantly different from that for identical atoms, equation (j^Hj) . For 
example, if the atoms are close together, 712 ~ 1, and then equation pTjl reduces to 



2N+1 ^ ' 

In this limit the concurrence is always positive, increases with A^ and approaches unity at 
a large A^. This is in contrast to the case of identical atoms where values of concurrence are 
below 0.25 even for 712 = 1 and approach zero for large A^. This behaviour can be easily 
explained by the fact that in the case of nonidentical atoms and 712 = 1 the population 
stored in the symmetric and antisymmetric states, pss + Paa, is equal to zero. At the same 
time, pu tends to unity as A^ increases giving the maximum concurrence Ci = 1. 

In real situations, we have 712 < 1, and for large A^ the concurrence Ci, given by (j2Z|), 
goes to zero similarly to the case of identical atoms. The maximum of Ci for nonidentical 
atoms, nonetheless, is much more pronounced than that for identical atoms. The real 
scale of large number of photons is in this case given by A^(A^ + 1)(1 — 7^2) ^ 1 rather 
than by A^(A^ + 1)(1 + 7^2) ^ 1 as it is the case for identical atoms. 

The dependence of the concurrence C = max(0,Ci) with Ci given by formula ()27j) on 
the interatomic distance ri2 and the mean number of photons of the squeezed field A^ is 
shown in Fig. |21 The dependence on the interatomic distance is similar to that seen for 
identical atoms with the revival of concurrence for ri2 ~ 3A/4 and not too large A^. For 
7i2 7^ 1 the dependence on A^ is also similar to that for identical atoms, except that the 
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Figure 2. Concurrence C for two nonidentical atoms as a function of the interatomic 
distance ryijX and the mean number of photons iV for |M| = y/N{N + 1) 

maximum values of concurrence for given are much higher than for identical atoms. 
Comparing the solutions for identical and nonidentical atoms indicates that one reason for 
higher values of concurrence for nonidentical as compared to identical atoms interacting 
with the squeezed vacuum is the fact that for nonidentical atoms less population remains 
in the lower block of the density matrix (jTHI) represented by states |3) and |4) (or \a) 
and |s)) as the atoms become different, i.e., A = {002 — uJi)/2 becomes large. To confirm 
this fact we plot in Fig. El the concurrence as well as the populations, p^s + Paa and 
Pgg + Pee, which are stored in the two blocks of the density matrix (fTH|l . as a function of A. 
Since for identical atoms, according to (0), considerable amount of population remains 
in the antisymmetric state in contrast to the solutions ()12j) for nonidentical atoms, it is 
clear from Fig. El that as the transition frequencies of the two atoms become more and 
more different the population of the antisymmetric state goes down reducing the total 
population pss + paa of the lower block and increasing the total population pgg + pee of 
the upper block of ()16p. which means higher values of concurrence. 

Another physical explanation of the origin of the better entanglement for non- 
identical atoms is provided by the observation that the stationary state of non-identical 
atoms, for small for which concurrence is maximal, is close to a pure state, whilst the 
stationary state of identical atoms is already far from a pure state. This is illustrated in 
Fig. m where we plot the purity measure P = Tr(p^) = + p^^ + p^^ + p^^ + |p„P/2 as 
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Figure 3. Concurrence C (solid line), populations pss + Paa (dashed line), and 
populations pgg + p^e (dashed-dotted line) for two atoms as a function of A for 
ri2/A = 0.05 and = 0.1 for |M| = sjN{N + 1) 




0.4' ■ ■ ■ ' 

0.5 1 1.5 2 

N 

Figure 4. Purity P = Tr(p^) as a function of the mean number of photons N for 
f'12/A = 0.05 and \M\ = y^N{N + 1): identical atoms (solid line), nonidentical atoms 
(dashed line) 



a function of N for the steady-state of two identical as well as nonidentical atoms. It is 
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seen that in both cases the purity decreases as the number of photons increases, but in 
case of identical atoms the purity goes down much faster. 




Figure 5. Fidelity F+ as a function of the mean number of photons A'' for ri2/A = 0.05 
and |M| = y/N{N + 1): identical atoms (solid line), nonidentical atoms (dashed line) 



It should be emphasised, however, that the main source of entanglement in the system 
are the nonclassical two-photon correlations that create two-photon coherences between 
the states \g) and |e). The two-photon coherences are nonzero only when the squeezing 
parameter \M\ is nonzero. In fact to have entanglement in the system the squeezed 
field must represent quantum correlations with \M\ > N. There is, moreover, one more 
necessary condition to have nonzero p„, which is a nonzero value of the collective damping 
parameter 712. The two photon coherences cause the system to decay into entangled 
states involving the ground state l^f) and the upper state |e) without any involvement of 
the entangled states \s) and \a). Unfortunately, the spontaneous emission from the state 
|e) redistributes some of the atomic population over the states |s) and |a) limiting in this 
way the degree of entanglement. 

Since the two-photon coherences create superposition (entangled) states involving 
the states \g) and |e), one can ask a question: How close is the entangled stationary state 
of the system to one of the maximally entangled Bell states 

l<^>^)= -^{\9) + \e)), 



1$ 



V2 



|e)). 



(29) 
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To answer this question we calculate the fidelities 

F+ = ($+|p|$ + ) = ^ {pgg + Pee + Pu) , 

F_ = ($-|p|<f-) = \ {pgg + Pee - Pu) ■ (30) 

The fidelities depend on whether pu is positive or negative. For small interatomic distances 
7i2 is positive, and then the coherence p„ is positive. Hence, the fidelity becomes 
large while the fidelity -F_ is small. Thus, we can conclude that the stationary state of 
the system is close to the maximally entangled Bell state |$+). In Fig. El we plot the 
fidelity as a function of for identical as well as nonidentical atoms. Comparing the 
dependence on N of the fidelity -F+ and concurrence C gives us clear evidence that the 
entanglement in the system can be related to the Bell state |$+). The state of the system 
is of course mixed, but it is closer to the pure Bell state |$+) for nonidentical atoms than 
for identical atoms. The entanglement created by the two-photon correlations present in 
the squeezed light is limited by the population stored in the other states. As the number 
of photons increases, for 712 < 1, more and more population goes to the states \s) and |a), 
and eventually entanglement disappears for both identical as well as nonidentical atoms. 
There are optimal values of the mean number of photons for which the highest possible 
stationary entanglement can be obtained. 

4. Conclusions 

In this paper, we have studied analytically the entanglement creation in a system of two 
atoms interacting with a squeezed vacuum field. We have demonstrated that nonclassical 
two-photon correlations characteristic of the squeezed field can create a large steady-state 
entanglement in the system. 

In our approach we have used master equation to describe a system of two two-level 
atoms subjected to a squeezed vacuum field. The two atoms are coupled to each other via 
the vacuum field which leads to the collective damping and collective dipole-dipole type 
interaction between the atoms. We have assumed the two atoms to be separated by a 
distance ri2, so the collective parameters depend explicitly on this distance. Steady-state 
solutions for the atomic density matrix have been found for two cases: (i) identical atoms, 
(ii) nonidentical atoms. 

We have derived analytical expressions for concurrence which is used to quantify 
the amount of entanglement created in the system. Our results show that the necessary 
condition for entanglement are nonclassical two-photon correlations of the squeezed field. 
The entanglement also depends on the interatomic separation and the mean number of 
photons of the squeezed vacuum. The necessary condition for entanglement are quantum 
correlations of the squeezed field. There is no entanglement for classically correlated 
field. We have found that the degree of entanglement created in the system is a result of 
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competition between the coherent process of transferring two-photon coherences from the 
squeezed vacuum to the atomic system and the incoherent process of spontaneous emission 
redistributing atomic population over the states not involved in the former process. In 
particular, we have shown that there is an optimum value of the mean number of photons 
for which the concurrence takes its maximum, and it happens for small number of photons. 
This is important from the point of view of practical applications. Moreover, we have also 
found that the degree of entanglement obtainable in this way is much higher when the two 
atoms are not identical. We have discussed in detail physical reasons for such behaviour 
of the two-atom system. 
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